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We study the asymptotic behavior of the eigenvalue distri- 
bution of the corner transfer matrix (CTM) and the density 
matrix (DM) in the density-matrix renormalization group. 
We utilize the relationship DM = CTM 4 which holds for non- 
critical systems in the thermodynamic limit. We derive the 
exact and universal asymptotic form of the DM eigenvalue 
distribution for a class of integrable models in the massive 
regime. For non-integrable models, the universal asymptotic 
form is also verified by numerical renormalization group cal- 
culations. 

PACS numbers: 05.50.+q, 05.20.-y, 75.10.Jm, ll.30.-j 



The density matrix renormalization group (DMRG) in- 
vented by S.R. White is one of the most important 
numerical methods developed recently. Due to the re- 
markable success, the method has now become one of 
the standard methods for studying one-dimensional (ID) 
quantum systems and two-dimensional (2D) classi- 
cal statistical systems. j| In spite of the success, lit- 
tle has been understood as for the foundation of the 
method. Studies clarifying the "origin" of efficiency of 
the method are important because they lead to vari- 
ous (including higher-dimensional 0) extension of the 
method. One example is the work of Ref. |5j, where a 
relationship between the density matrix (DM) and Bax- 
ter's corner transfer matrix (CTM) [Q^is pointed out 
and a new algorithm (CTMRG) is devised. Another ex- 
ample is the work of Ref. pfj| , where it is pointed out 
that the DMRG (at its thermodynamic limit) is a vari- 
ational method using the matrix-product-ansatz (MPA) 
wavefuntion as a trial wavefunction. This leads to a di- 
rect variational method which does not need the DM [jlOj , 
and the product- wavefunction RG (PWFRG) which fully 
utilizes the MPA form of the DMRG-fixed-point wave- 
function, n 

The central object in the DMRG is the DM which is 
made from the groundstate wavefunction (resp. maximal 
eigenvalue wavefunction) of quantum Hamiltonian (resp. 
transfer matrix) by tracing out information of one half of 
the system. Keeping up to a cut-off-eigenvalue eigenstate 
of the DM, we have a truncated basis set consisting of a 
finite number (conventionally denoted by "m") of bases, 
to describe the remaining half of the system. 

Since the accuracy of the DMRG is determined by the 
cut-off eigenvalue, it is crucially important to investigate 



the eigenvalue spectrum {w m }, in particular, its asymp- 
totic (m — > oo) behavior which has not been known pre- 
cisely. In this Letter, we present the exact asymptotic 
form of the DM eigenvalue distribution for a class of 
(non-critical) integrable models, and further, make the 
first systematic study for non-integrable systems employ- 
ing the CTMRG and the PWFRG by which wc can effi- 
ciently obtain the "fixed-point" (thermodynamic limit of 
the system) of the DMRG. 

Let us first discuss the integrable cases. In these cases, 
ID quantum problems are equivalent to 2D classical sta- 
tistical problems: the Hamiltonian of the former can be 
derived by a log-derivative of the transfer matrix of the 
latter, and the ground-state wavefunction of the former is 
identical to the maximal-eigenvalue eigenfunction ^/ max 
of the latter. M Hence, we discuss only 2D classical cases 
in the below. 

As has been pointed out by Baxter, js| the wavefunc- 
tion (WF) f mas is interpreted as a product of two CTMs, 
in the thermodynamic limit. Since the DM is just a 
square f^ lffl (with f ma x being regarded as a "wavefun- 
tion matrix"), this interpretation leads to a relationship 
H between the DM, the WF and the CTM for 2D classi- 
cal systems (at least non-critical case where the boundary 
effect is negligible), which is symbolically written as 



WF = (CTM) 2 , 
DM = (CTM) 4 . 



(1) 



For integrable models, diagonal form of the CTM is easily 
known, from which we can obtain, for example, the exact 
one-point function (spontaneous magnetization, etc), j^] 
Due to the relationship ([!]), the diagonal form is also 
useful to obtain the exact eigenvalue spectrum of the DM. 

We discuss the simplest case where the diagonal form 
of the CTM is given by a single infinite tensor product; |jj 
due to ([!]) , diagonal form of the DM has the same infinite- 
tensor-product form with redefined parameter. For def- 
initeness let us consider two cases (Type I and Type II) 
where the exact diagonal form of the DM is given by 
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(diag) _ 







(2) 



with c„ = n for Type I models (e.g, transverse-field Ising 
chain, 6-vertex model, eight-vertex model, XXZ-chain 
and XYZ-chain) and c n — 2n — 1 for Type II models 
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(e.g., the square-lattice Ising model in the conventional 
(not eight- vertex) representation). B The parameter z 
(0 < z < 1) represents "degree of non-criticality" (i.e., 
z — > 1, on approaching the critical point), and how it 
relates to "physical" parameters depends on the model. 
Note that the DM (||) is unnormalized. It is "normal- 
ized" in such a way that its maximal eigenvalue ujq is 
unity; we should divide it by Tr/3( dlag ) for conventional 
normalization. 

Due to the tensor-product structure (||), each eigen- 
value of p( dlag ) has the form z n with n (> 0) being an 
integer. Further, each eigenvalue z n may have degener- 
acy p(n). To study the degeneracy structure of the DM, 
it is convenient to consider Trp dlag : 



Tr/ iag = \{(l + z^) = Y J P(n)z n , 



(3) 



n=l 



ri=0 



where the degeneracy p(n) is precisely the coefficient of 
z n in the infinite series. We should note that, taking the 
degeneracy into account, the number of retained bases m 
in the DMRG should be 



m = m(n) — J^p(fc), 

k=0 



(4) 



which means that the cut-off eigenvalue of (unnormal- 
ized) DM is z n and that we should retain all the degen- 
erate bases belonging to this cut-off eigenvalue. 

Our problem is to obtain the large-n behavior of m = 
m{n). For this purpose, we should know the asymptotic 
behavior of p(n). The partition theory of integers, which 
has been played an important role for studies of inte- 
grable IRF (interaction-round- a- face) models, is helpful 
again. | pT| By r(n) we denote the number of partitions of 
a positive integer n under a restriction "r" Consider the 
generating function f(q) associated with the restricted 
partition problem. It has been known Jj"l| that for a 
wide class of partition problems, f(q) can also be given 
in an infinite-product form: 



n=0 n=0 



(5) 



where each a n is a non-negative real number. For the 
Type I case, we have a n = 1 for n odd, and a n = 
otherwise. 

The asymptotics of the generating function of the form 
(||) is calculated by the saddle point method; r(n) for 
ji > 1 is then given by Meinardus's theorem (cited in 
Ref. [0, page 89): 



~(n) = An K exp(Bn a ^ 1+aS> ) + (less dominant terms) 



oo 
n=l 



and k is given by 



K = 



D(0) - 1 - a/2 
1+a 



(7) 



(8) 



Explicit forms of A and B which we have omitted in the 
above are also given by the Meinardus's theorem. For 
the Type I models, we have a = 1 and k = —3/4: 



p(n) = const n 3 / 4 e?q>(B^/n), 



(9) 



where B — ir/y/3. |Tl[] For the Type II models, a related 
theorem (Ref. [0, page 99-100, example 10 and 11) as- 
sures the same asymptotic form (||) with B — it/ \/6. It 
is also possible to relate Type II models with the Meinar- 
dus's theorem (Chap.l and 6 in ref. 0). We thus have 
derived the exact asymptotic form of p(n) . 

Using (|4|) and changing the summation into the inte- 
gration, we finally obtain 



,-1/4 



exp(By/n), 



(10) 



for the Type I and II models. How well the DMRG 
calculation for the 5 = 1/2 XXZ chain reproduces the 
asymptotic behavior ([Toj ) is demonstrated in Fig.l. In 
the actual calculation, we have employed the quantum 
version of the PWFRG, |p^-p^ by which we obtain the 
fixed-point wavefunction of the DMRG efficiently. 

We give a comment on the universality of the asymp- 
totic form ( |l0| ) among the integrable systems. In the case 
where {a n } forms a periodic series or the model itself 
admits a direct partition-theoretic interpretation, [jl5 16 
the exp(By / n)-behavior is universal (Ref. pj . Chapter 
6, examples 1-16). The exponent k may, however, have 
possibility of model-dependence (due to D(0)), modify- 
ing the prefactor n -1 / 4 in ([l0|). 

Let us now proceed to non-integrable cases, where the 
exact diagonal form of the CTM or the DM is not known. 
The DM eigenvalue is no longer given by z lnteger with 
single parameter z, or equivalently, log(DM eigenvalue) 
has not equal-spacing distribution. Both the integer n 
characterizing the DM eigenvalue, and the quantity p(n) 
which is essential in the integrable cases lose meaning. 
Our first task is, then, to translate the result of inte- 
grable cases into the one which has meaning also for non- 
integrable cases. 

Writing the m-th DM eigenvalue (including degener- 
acy) as w m , we have n = logw m /logz in the integrable 
case. Substituting n = logw m /logz into (|l0|), we have 



(6) 



logz 



-1/4 



exp 



nog to m 
logz 



(11) 



where a is the real part of the pole of the Dirichlet series, or equivalently, 
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log 



logz 



1/4' 



= B\ 



< \QgLO m 

\ogz 



(12) 



From (|l2|), we obtain the leading asymptotic form 

u> m ~ exp[— const, (logm) 2 ], (13) 

where const. = |logz|/£> 2 for the intcgrablc cases. 
Clearly, expressions (|Ii])-(|l]) do not contain the param- 
eter n which is specific to the integrable models. 

There arises an intriguing conjecture: the asymptotic 
forms (|ll|)-(|l3|) would also apply to non-integrable sys- 
tems with B and z being suitably redefined. In the 
"neighborhood" of an integrable model with small non- 
integrable perturbations added, we may well expect this 
conjecture to be true: In spite of the non-integrable per- 
turbations, the "stairway structure" (or degeneracy) in 
the DM eigenvalue spectrum still remains in somewhat 
smeared-out way, leaving the "envelope" of the Lo m -m 
curve essentially unchanged. As a check of the uni- 
versality for the nearly-integrable cases, we made the 
CTMRG calculations for two systems: the square-lattice 
Ising model at the critical temperature in finite exter- 
nal field and the 3-states Potts model slightly below the 
critical temperature (see Fig. 2 and Fig. 3). We see clear 
agreements between the CTMRG calculations and the 
"universal asymptotic form" . 

As a test of the universality of (|Il|)-(|l3|) for systems far 
from the integrability, we take the S = 1 antiferromag- 
netic Heisenberg spin chain. For calculation of the DM 
eigenvalue spectrum, we employ the quantum version of 
the PWFRG. pip) The results are given in Fig.4, which 



support the universal asymptotic form. 

We have made similar calculation for the 5 = 1 
bilinear-biquadratic spin chain at — —0.5 with the 
Hamiltonian H = Y^,^i ' Si+i + PY^i^i ' Si+i) 2 , whose 
result (not shown in this Letter) also supports the uni- 
versality of the asymptotic form. 

To summarize, we have discussed the asymptotic dis- 
tribution of the density-matrix (DM) eigenvalues for 
non-critical systems (one-dimensional quantum and two- 
dimensional classical), which controls the accuracy of 
the density-matrix renormalization group. Based on the 
equivalence between the DM and the corner transfer ma- 
trix (CTM), we derived the exact asymptotic form of 
the DM eigenvalue distribution for a class of integrable 
models. The resulting expression has been rewritten in a 
"universal" form which does not contain quantities spe- 
cific to integrable models. Numerical-renormalization- 
group calculations using the CTMRG and the product- 
wavefunction RG have been performed for non-integrable 
models, which shows that the non-integrable models 
actually have the same asymptotic form of the DM- 
eigenvalue distribution, in strong support of the univer- 
sality of the asymptotic form. 

There remains many important problems left for fu- 
ture studies. A more "physical" explanation to justify 



the universal asymptotic form is desired. How univer- 
sal the obtained asymptotic form, itself remains to be a 
question to be answered; there may well be different "uni- 
versal classes" of the DM. In fact, the valence-bond-solid 
(VBS) models [|l7|, have only finite-dimensional DMs 
which sharply contrast to the ones studied in this Letter. 
Relation between the DM-eigenvalue distribution and the 
finite-m (m: number of retained bases) behavior of physi- 
cal (observable) quantities has not been known, although 
there have been a few works discussing the "finite-m scal- 
ing" (2D classical @, transverse-field XXZ chain 
Behavior of the DM for critical system is also an im- 
portant subject of study. J2(],|2l| Our study made in the 
present Letter may be a first step for clarification of these 
problems. J2^] 
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FIG. 1. PWFRG calculation of density-matrix eigenvalues 
{u; m } for S = 1/2 antiferromagnetic XXZ chain and compari- 
son with the exact spectrum. We take the exchange coupling 
constants to be | J x \ — \ J y \ = 1 and \J Z \ = A = cosh(l). Num- 
ber of retained bases in the PWFRG calculation is m — 207. 



FIG. 2. CTMRG calculation (m = 200) of the den- 
sity-matrix eigenvalues {u) m } for the square-lattice Ising 
model at a critical temperatur T c in a small external field 
H. We have also drawn a line corresponding to the universal 
asymptotic form. 



FIG. 3. CTMRG calculation (m = 242) of the den- 
sity-matrix eigenvalues {tj m } for the 3-state Potts model 
slightly below the critical temperature. We have also drawn 
a line corresponding to the universal asymptotic form. 



FIG. 4. PWFRG calculation (m = 700) of density-matrix 
eigenvalues {um} for S — 1 antiferromagnetic Heisenberg 
chain. We have also drawn a line corresponding to the uni- 
versal asymptotic form. 
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